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Abstract 

The dynamics of a heteropolymer chain in solution is studied in the limit 
of long chain length. Using functional integral representation we derive an 
effective equation of motion, in which the heterogeneity of the chain man- 
ifests itself as a time-dependent excluded volume effect. At the mean field 
level, the heteropolymer chain is therefore dynamically equivalent to a ho- 
mopolymer chain with both time-independent and time-dependent excluded 
volume effects. The perturbed relaxation spectrum is also calculated. We 
find that heterogeneity also renormalizes the relaxation spectrum. However, 
we find, to the lowest order in heterogeneity, that the relaxation spectrum 
does not exhibit any dynamic freezing, at the point when static (equilib- 
rium) "freezing" transition occurs in heteropolymer. Namely, the breaking of 
fluctuation-dissipation theorem (FDT) proposed for spin glass dynamics does 
not have dynamic effect in heteropolymer, as far as relaxation spectrum is 
concerned. The implication of this result is discussed. 
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I. INTRODUCTION 



Although systems without quenched randomness usually provide useful information and 
serve as the first step of our understanding of condensed phase, there are situations in which 
quenched randomness brings entirely new phenomena. Flux line pinning in type II super- 
conductors, critical to the practical application of superconductors in high field magnets, is 
a manifestation of the embedded quenched impurities M. In the mixture of 3 He and 4 He, 
quenched randomness introduced by porous media (e.g. aerogels) shifts the tricritical point 
and extends the range of superfluid phase M. The interesting properties (e.g. the onset of 
remanence effects below a freezing temperature) of dilute magnets, modeled by "spin glass" 
||, are consequences of the magnetic impurities. These are just a few examples of the new 
phenomena brought by quenched randomness. Study of these phenomena constitutes an 
important part of modern condensed phase physics. 

Theoretical progress in the study of these systems, especially spin glass systems, also has 
brought impetus to other, somewhat distant, fields. For example, it is shown that there is 
an interesting analogy between some models developed for spin glasses and the statistical 
properties of polymers with quenched randomness (realized as heterogeneity of the polymer) 
This analogy has attracted much attention in recent years due to its possible biophysical 
applications. As much work has been focused on the static (thermodynamic) consequences 
of this analogy, we shall discuss in this paper its dynamic consequence. 

One of the general methods used to study systems with quenched randomness is the 
so-called replica method HHH. Because the randomness is quenched, thermal equilibrium is 
established without thermally equilibrating with the randomness. Therefore, to calculate 
the thermodynamic quantities for a system with quenched randomness one has to calculate 
the average of the free energy over the distribution of the randomness, 

(-k B T\nZ). (1.1) 

The mathematical obstacle lies in the difficulty in calculating the average of a logarithmic 
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function depending on the randomness. Replica method circumvents this obstacle by using 
the identity 

Zi n — 1 

lnZ = lim (1.2) 

n-»0 n 

Replica method has been very useful and successful, though not free from controversy. 
With its success and popularity in spin glasses, it is not surprising that the study of het- 
eropolymers along the line of the above-mentioned analogy often, if not exclusively, uses this 
method to draw thermodynamic conclusions. 

Besides the replica method, there exist other useful theoretical methods to study sys- 
tems with quenched disorder ||. Since we are interested in the dynamics of heteropolymers, 
among these methods the Langevin dynamic method |||6j is the most relevant one. This 
method does not use the unphysical replica limit n — » to calculate thermodynamic quan- 
tities. In Ref. it is shown that the functional integral formalism developed by Martin, 
Siggia, and Rose for classical statistical dynamics [0] allows average over the quenched ran- 
domness without using the replica method. It is also shown that this formalism is capable 
of giving dynamic information as well as statics. Therefore, although it generally involves 
more complicated techniques, it is worthwhile to study systems with quenched randomness 
using this method. In spin glasses there have been some studies using this approach [|]|§J1. 
On the side of heteropolymers, however, to the best of our knowledge, a dynamic study 
similar to these studies is still lacking. While the static analogy between spin glasses and 
heteropolymers has been very useful, a similar dynamic study is necessary if one wants to 
use this analogy to address dynamic issues. 

The purpose of this paper is to study the dynamics of single heteropolymer chain, i.e. 
dilute heteropolymer solution. It is well known that the dynamic study of polymer solu- 
tions, especially dilute solutions, must take into account the hydrodynamic effect in order 
to make the prediction experimentally relevant [pl3|-|T2"|j . Besides hydrodynamic effect, the 



heterogeneity of the polymer chain presumably introduces further effect. One expects that 
the dynamics of the heteropolymer chain will depend on these two effects as well as the 



usual excluded volume effect. Here we investigate the importance of the heterogeneity to 
the chain dynamics. 



The dynamics of dilute polymer solution has been extensively studied [ 10,12]. Dynamic 



quantities such as diffusion constant ||13|| , time- correlation functions [Tl], intrinsic viscosity 



initial decay rate of the dynamical scattering and relaxation spectrum [P7|] have 



been calculated. Among these dynamic quantities we are most interested in the relaxation 
spectrum because, in the context of Ref. [Q, i.e. protein dynamics, the relaxation spectrum is 
the most relevant quantity with potential application to the protein folding problem. In spite 
of the close analogy between models of spin glasses and heteropolymers at static level 0, we 
find that this analogy is not as complete at dynamic level. While there have been attempts to 
study heteropolymer dynamics at the phenomenological level of the hypothetical dynamics 
of the Random Energy Model (REM) |18[ , the present microscopic analysis suggests that 
such analogy may not be justified when dynamic issues are tackled. 

The structure of this paper is as follows. In Sec. [TT] we define the model Hamiltonian 
and the dynamic equation of the system. We discuss the physics of the model Hamiltonian 
and the dynamic equation. Assumptions made are also addressed. In Sec. |TTT| we derive the 
mean field equation of motion for the chain in the limit of long chain length. This effective 
equation of motion can be further simplified when the characteristic time scale of the solvent 
is much shorter than the characteristic time scale of the chain. An effective Hamiltonian 
results from this effective equation. It is shown that the effective Hamiltonian contains a 
memory term, and, hence, a time-dependent excluded volume interaction. In Sec. [IV A 



we 



calculate the perturbed relaxation spectrum for the chain, to the lowest order in excluded 
volume interaction, hydrodynamic interaction, and heterogeneity. We find that, as long as 
the model Hamiltonian and dynamic equation are valid, the relaxation spectrum does not 
depend on the final equilibrium state of the chain. In other words, in spite of the static 
analogy between spin glass and heteropolymer mentioned above, the possible breaking of 
FDT, due to the degeneracy of the ground states, found in the mean field dynamics of spin 
glass UK does not have corresponding dynamical effect, at least to order of our calculation, 



in the mean field dynamics of heteropolymer. In Sec. |V| we discuss the analogy of spin glass 
and heteropolymer, using the well-known correspondence between magnetic system and 
self- avoiding random walk. We will see that, when the kinetic equations are compared, it is 
clear that the analogy is not very helpful, because the dynamic variables used are completely 
different and may not be relevant to each other. In Sec. [VI] the physical situations to which 
our model may be applicable are discussed. We also give a brief discussion of approximations 
involved and the limitation of the results. In App. [A| we provide an alternative derivation of 
the conventional (Kirkwood) approach used in the dynamic study of dilute polymer solution. 



App. H gives the detailed derivation of the mean field approximation sketched in Sec. |TI 
In App. y we detail the calculation leading to the result in Sec. [IV A 



II. THE MODEL 



A. Hamiltonian 



Consider a heteropolymer chain with conformation specified by the coordinates of N 
monomers, ro, • • • , ?n-i- The chain can be described by the following Hamiltonian, 

H = H + H R , (2.1a) 

N-l i 

H = £ -(r m - ri ) 2 , (2.1b) 

i=o z 

Hr^Wb^U^-v^, (2.1c) 
where B^ is a Gaussian random number centered at Bq with variance B 2 (B > 0), 
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V2^B~ 2 



(Bij - Bq) 
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P (^) = ^™ ex P ~ 2B 2 ' (2 ' 2) 



In Eqs. (|2.1| ) Hq is the Hamiltonian used in the bead-spring chain model (Rouse model) 
where, for convenience, all the spring constants have been taken to be unity (see below). 



The chain connectivity of the polymer is described by this term. Hr simulates the excluded 



volume effect. The excluded volume interaction U (r* — r^) between monomer i and monomer 
j is Gaussian-modulated by Bij. This modulation confers the chain heterogeneity. When 
B —>■ the excluded volume interactions are "monochromatically" modulated by Bq and 
the chain becomes homogeneous. At this extreme, the Hamiltonian, Eq. ( |2.1| ), reduces to 
the usual Edwards Hamiltonian JT3] , 



N-l 1 D 

i=0 i^j 

On the other hand, a nonvanishing means that there exist various monomer-monomer 
interactions, and the chain is essentially heterogeneous. Thus B is a measure of chain 
heterogeneity and will be called heterogeneity parameter. Bq is the usual excluded volume 
parameter measuring the strength of excluded volume effect. Presumably, B depends only 
on the heterogeneity (dispersity) of polymer, while B is generally a function of temperature 
and solvent quality. 

The excluded volume parameter B will be taken to be positive throughout this paper. 
A negative Bq gives attractive monomer-monomer interactions, which requires inclusion of 
repulsive three-body interactions in order to render the theory stable |1^,[2(J. To simplify 
the matter, we will confine ourselves to positive B , i.e. a heteropolymer chain in good 
solvent. 

The conventional use of the excluded volume interaction U(r) is the hard-core repulsion 
5(r) |T2]JT9|j20|] . To derive the effective equation of motion at the mean field level in a more 



general form, we will not adopt this convention for U(r) until Sec. [IV A| , where the perturbed 
relaxation spectrum is calculated. Before Sec. [IV A| the only assumption on U(r) will be 
translational invariance: U(r, r') = U(r — r') = U(r f — r). 

Obviously, the Hamiltonian H ignores the possible variation of spring constants. Pre- 
sumably, this variation exists when the chain is heterogeneous. Since our concern here is the 
effect of random excluded volume interactions, we do not consider this variation in spring 
constants. With this simplification, by a properly chosen length scale, we can always set all 
spring constants to be unity, as is done in (|2.1b|) . 
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The statics of Hamiltonian (|2.1|) has been studied in Refs. |4]], within the context of 
biopolymer and protein folding. These works share a number of techniques with spin-glass 
models. In the past two decades there has been much progress in the theories of spin-glass, 
in both statics and dynamics || . However, there has never been any systematic microscopic 
study on the dynamics of Hamiltonian ( |2.1[ ). 



B. Dynamics 

Starting with the Hamiltonian (|2.1|) , for dynamic study, we consider the coupled Langevin 
equations proposed by Oono and Freed [pi]] : 

d 1 dH 

—Ti(t) = --——+g oU (T i (t),t)+0 0i (t), (2.4a) 
dt Co OTi(t) 

(0 0i (t)) = 0, (0 O i(*)MO> = T 5{ 1 ~ *')M> ( 2 - 4b ) 

a _ N-l air 1 

£u(r, t) = u Au(r, t) - °* £ -^(r - r 4 (t)) - - Vp + ^ (r, t), (2.4c) 

(£o(r, t)> = 0, (£ (r, t)^ (r', 0> = -^?A5(r - r')5(t - t')l, (2.4d) 

Po 

where 1 is a unit tensor. In these equations the kinetic coefficient (q 1 and kinematic 
viscosity v$ set the time scales for polymer and solvent molecules, respectively, go is a 
measure of the coupling strength between monomers and solvent molecules. (We also ignore 
the possible variation of these couplings due to the chain heterogeneity.) Langevin noise 
0Qi(t) gives monomer % a random velocity at time t. Noise £ ( r ' t) ^ s ^ ne random acceleration 
of the solvent velocity field at position r and time t. Dynamic viscosity 7] is, as usual, 
related to the solvent density p by r] = po^o- The hydrostatic pressure p ensures the 
incompressibility of the solvent ||22|| . This condition enables us to consider only the transverse 
component of solvent velocity field. Since no confusion will be arisen, hereafter, we will 
denote the transverse component of solvent velocity by the same symbol u. 



Eq. ( [2.4a]) is a relaxation dynamic equation, with mode-coupling term goii, for the chain. 



Eq. ( |2.4c| ) is the Stokes approximation of the Navier-Stokes equation for solvent velocity 
field u, augmented by the forces exerted by monomers [^] . We will assume that the spatial 
extension of the system is infinite so that we can ignore the boundary condition. The problem 
related to the non-applicability of the Navier-Stokes equation in an unbounded space will 
be also ignored |]23|^4]| . Eqs. (|2.4j) are valid as long as the Reynolds number is small. 

In the conventional theory of Brownian motion the effect of surrounding solvent molecules 
is taken into account by the Langevin noises. In polymer dynamics the commonly used ap- 
proach (Kirkwood's approach) is essentially based on this viewpoint [|T~0| P~2|] . This approach 
treats each monomer as a point source of frictional force and lumps up the effect of sol- 
vent molecules in the Langevin noises. However, as argued by Oono |T^ , |25[| , this picture 
is not self-consistent because the monomers are not legitimate Brownian particles, even if 
the polymer as a whole can be treated like a Brownian object. In contrast to Kirkwood's 
approach, the coupled Langevin equations Eqs. ( |2.4| ) proposed by Oono and Freed explicitly 
introduces solvent velocity field. In this approach the Langevin noises 6q an d £o have to be 
re-interpreted as coming from the coarse-graining procedures used to derive these equations, 
separately, from the more fundamental, microscopic, equations, and the coupling between 
polymer and solvent is introduced only after the respective coarse-graining procedures have 
been performed. A proper term for these equations is "kinetic equation" ||26|| . Therefore, 
we will simply call Eqs. ( |2.4|) Oono-Freed kinetic equation. 

In the Oono-Freed kinetic equation the energy scale has been chosen as k^T = 1. The 
correlation of Langevin noises chosen here, Eqs. ( |2.4b| ) and ( |2.4d|) guarantees that the system 



will eventually approach to the equilibrium state determined by the statics of Hamiltonian 
Note that the second term on the right side of Eq. ( |2.4c| ) is the frictional forces exerted by 



the monomers, which are point sources of frictional force, same as in Kirkwood's approach. 
In fact, there is a close relationship between these two approaches. It has been shown that, 
to the lowest nontrivial order (i.e. to the lowest order in the excluded-volume parameter Bq 
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and hydro dynamic coupling parameter go) and within Markovian approximation for solvent 
velocity field, when the solvent velocity field is projected out, the Fokker-Planck equation 
for the kinetic equation Eqs. ( |2.4j ) reduces to the conventional Kirkwood diffusion equation 
I2| , [l6|| . A simplified proof of this reduction is given in App. [A]. These equations are, in this 



sense, more fundamental than the commonly used Kirkwood diffusion equation, upon which 
most dynamic study of dilute polymer solution are built. Therefore, a full dynamic study 
should start from Eqs. (|2.4| ). 

Using the same spirit of the proof in App. [A], we first project out the velocity field from 
the kinetic equation. This is done by formally solving Eq. (|2.4c ), 



u(r,t) = u (r,t) + u R (r,t) 



Uo(r,f)= / e lk - r u (k,t) 



JV-l 



dH, 



Jk J-oo L Po f^Q OTiiif) 

u R (r,t)= I e 4k - r u (k,t) 



o_ e -ik. ri (t') 



r ft r n N ~ 1 

ik J — oo L On ._n 



9o U dH R 



-ikTi(t') 



(2.5a) 



(2.5b) 



(2.5c) 



Po t^o 5r *( f ) 

In Eqs. ( |2.5b|) and ( |2.5c|) P T (k) is the transverse (along k = \t/k direction) projection 
operator, 



P T (k) = 1 - kk 



and 



d d k 



(2.6) 



(2tt) c 



(2.7) 



(d is the spatial dimensionality). The Fourier-transformed random acceleration, 

is Gaussian and is correlated according to 

(^(k, t)^(k', £')> = P T (k)^A; 2 5(k + k')5(t - f ) 



(2.8) 



(2.9) 



where _L denotes the transverse component, £ (k, t) = P T (k)-£(k, t). 

Substituting solution Eq. (|2.5| ) into Eq. ( |2.4a|) , the equation of motion for the chain 
becomes 

1 OH 1 OH 

-Mt) = --—°+g u o {T i (t),t)+e 0i {t) - -—JL+g U R (Ti(t),t), (2.10) 

at ^OYi[t) ^QOVi[t) 



Eq. ( |2.10| ) is the starting equation of our dynamic study. Note that this equation has two 



stochatic terms: g u (ri(t),t) which contains £ (k, f), and Ooi(t). In other words, the new 
Langevin noise is 

e 0i (t)+9o I f dt'e-»° k2 ^P T (kH (k,t') (2.H) 

J\l J — oo 

Obviously, this new Langevin noise is no longer 5-correlated in time. 

Eq. (|2.10|) is much more complicated than its spin-glass counterpart |||||J, because 
of the presence of the excluded volume interactions and the hydrodynamic coupling. In 
polymer dynamics, quenched randomness has been considered (for example, the mean field 
dynamics of directed polymer in random meida studied by Vilgis ||271|), but the randomness 
is from external random potential rather than chemical dispersity. Also, equations similar 
to Eq. ( 2.10] ) in which hydrodynamic effect is included have been studied in the context 



of homopolymers P^|-[T^j2~I|j2^|PB|j2^|| and two dimensional membranes ||30|| . These stud 



ies, however, consider simpler situations in which either quenched randomness or excluded 
volume interaction is absent. Eq. (|2.10|) considers all three effects: quenched randomness, 



hydordynamic coupling, and excluded volume effect. 

III. MEAN FIELD APPROXIMATION AND THE EFFECTIVE EQUATION OF 

MOTION 

A. Functional- integral representation 



The equation of motion, Eq. ( |2.1U| ), is essentially a (nonlinear) stochastic differential 
equation with colored noise (Eq. ( [2.11D ) [31| . A standard and convenient method to study 
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stochastic differential equations is to use functional integral formalism (Martin-Siggia-Rose 



(MSR) formalism) |7|J3T|-|33| . The main idea is to write down the probability functional for 
the Langevin noises and make change of variables from these noises to dynamical variables. 
The stochastic differential equation itself is treated as a constraint, limiting the evolution 
of probability path. By introducing auxiliary fields conjugated to the dynamical variables, 
this constraint, in the form of a 5 functional, can be written as a functional integral of these 
auxiliary fields. This leads to a probability functional in terms of the dynamical variables 
and their conjugated auxiliary fields. The auxiliary-field technique used in this formalism 
is similar to the technique used in the supersymmetry formalism of stochastic differential 
equations, in which the auxiliary fields are replaced by fermion fields [ |3~2| , j34}1 . 

Since the functional integral formalism for stochastic processes is already well docu- 
mented, we would not go into any detail of it. Useful references are Refs. |]7|, |32|j33[| . (The 



concise review in Ref. J35] is also helpful.) The MSR formalism has been widely used for 
dynamic study of condensed phase systems, for example, liquid-glass transition P5[l , tur- 
bulence ||36|| , spin glass dynamics || and so on. Nevertheless, in polymer physics it is not 
often used. This is probably because that, in spite of its elegance, no new results have 
been obtained through using this formalism |J7|,[3S|. The study in Refs. [JT5|jT^ ,P 



makes 



it clear that, for practical calculation, the MSR formalism does not achieve, at this point, 



better (higher order) results than the conventional approach |[T^j25|1 . The mathematical in- 
volvement makes higher order calculations insurmountable. However, for our purposes the 
MSR method is very useful as it allows to carry out quenched averaging in the most natural 
form (see below). 

Using the MSR formalism, the generating functional for Eq. ([2.10|) is presented as: 

Znhto = I {I?r J (t)}{I?f l (t)}J({r J })e i ° +L «, (3.1) 



N-l 



/oo ' ' * 
dt ™i(t) 
-°° i=0 

roc r \ QJ{ R 

Lr= dt22 i$i(t) 

J —OO ■ n 



d 1 dH 



g u R (ri(t),t) 



(3.2) 
(3.3) 
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Note that the Langevin noises 6q and £ are now in Lq, while the quenched heterogeneity, 
represented by B { ^ is in L R . Since , £ and are not coupled, we can take average of 
Zjie £ over these random quantities separately 

The Jacobian J({rj}) in Eq. ( |3.1| ) is associated with the change of variables from Langevin 
noises to dynamical variables (mentioned above). Because of the following two reasons, we 
will drop the Jacobian from our calculation from now on. Firstly, the specific form of 
this Jacobian depends on the time-discretization scheme used when the functional integral 
representation is written down f39|. (For Stratonovich scheme Refs. [fjOifldj work out the 
detailed form.) Although it is not uniquely defined, one can get correct result provided one 
consistently uses the same discretization scheme [f42[ . Secondly, and more importantly, the 
Jacobian depends on the dynamical variables {r^}, not on the auxiliary variables {?»}. As 
we can see from Eqs. ( |3.2j ) and ( |3.3|) the equation of motion, Eq. (|2.10| ), is always coupled 
to the auxiliary fields {?«}. This will still be true for the effective equation of motion that 
we seek for. Therefore, whatever form the Jacobian has, or, whether it also gets averaged 
when we perform averages over B^, 0q and £ is not important. We do not need it anyway. 

In ( |3.1| ) only the auxiliary fields {fj} are introduced. In principle, for the Oono-Freed 
kinetic equation fl2.4j), there should also be auxiliary fields {u}, as in Ref. (TBI.ini.pRpP 



(3.4) 



Since we are not interested in the solvent velocity fields, u, and have eliminated it from the 
kinetic equation, there is no need to introduce an auxiliary variable u for u. 
For computational convenience we define three operators, 

<i 

Po 

[a, symmetric tensor), 

1 r N ~ x r i 

a t (t) = - dt"£ iW> - t' ~ ~ + Tjitf^W ~ *) (3.5) 

1 J j=0 

(a scalar differential operator; e being an infinitesimal positive number), and 



r«(t,f) = ^Co / e -^ fc2 ( t -'')e ik -^w- r ^*')]p T (k) 

Po 



o 



dt 



ai {t)U{v t {t) - Tj (t)) + a^Uiv^t) - n(t)) 



(3.6) 
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(symmetric in % and j). Now we can rewrite Lr in a more compact form 

The average over B^ can be easily done and gives 
Z e ^ = j{Dv t {t)}{Di t {t)}e L »(e L *) B 

= [ {Vv i (t)}}{Vv i (t)}e Lo exp \* £ 0% + ^- y £O i 



(3.7) 



(3.8) 



B. Mean field approximation 



The generating functional obtained in Sec. [Ill A| can be used to derive an effective equa- 
tion of motion in the limit of long chain length N — > oo. The same method is used to derive 
an effective equation of motion for spin glass |||| and for directed polymers [|27|. Although 
the algebra involeve here is more complicated, the resulting effective equation of motion 
share the same characteristics as the effective equation of motion obtained in Ref. |||9]]. 

The detailed calculation in the long chain limit can be found in App. |B[ The result is 
an effective Lagrangian: 



B 2 

L +-7r I dtdt' dld2d3dA 



Co 2 



2Bn 



A(t, t', 1, 3))B(t, t', 2, 4) + (B(t, t', 2, 4))A(t, t', 1, 3) 



+2(C(t,t',l,3))D(t,t',2A) 



C/(l-2)[/(3-4) 



dtdt' dld2d3dA 



C(t,t',l,3))B(t,t',2,A) 



B(t,t',2,4))C(t,t',l,3) 



5(t-t')U(l-2), 



(3.9) 



where, for convenience, Ri, R 2 , R3, R4 are written as 1, 2, 3, 4 respectively. The quantities 
A, B, C, D are (see App. @) 

A(t, t\ R, R) = a^a^diviit) - R)%(t') - R) (3.10a) 

i 

B(t,t', R,R') =J2 5 ( r i(t) ~ R)S(ri(t') -K') (3.10b) 

i 

C(t, t', R, R) = ]T Oi(t)S(Ti(t) - R)S(ri(t') - R) (3.10c) 

i 

D(t, t', R, R) = ]T <kW)S(ri(t) - R)6(Ti(lf) - R) (3.10d) 
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Note that B is the dynamic version of the order parameter used in Ref. |4j]. In spin 
glass dynamics it is also found that the static Edwards- Anderson order parameter should be 
generalized to a time-dependent order parameter whose long time limit gives the static order 
parameter PJ4"3|. Similarly, the static order parameter used in Ref. can be defined as 
the long time (equilibrium) limit of this dynamic order parameter B. The physical meaning 
of the static order parameter is explained in Ref. {§]. We will explain later the physical 
meaning of the dynamical order parameter B. 

A little reflection tells us that, firstly, the terms {A)B and (C)B will not appear in the 
final effective equation of motion because they are related to the Jacobian when writing the 
functional integral representation for the effective equation of motion. As explained earlier, 
the Jacobian does not play an important role, as far as the effective equation of motion is 
concerned. Hence, we will also ignore these terms from now on. Secondly, the term (B)A 
is related to the stochastic part of the effective equation of motion because it contains two 
auxiliary fields outside of the self-consistent average bracket (). Finally, because of having 
only one auxiliary field outside of the self-consistent average bracket, the terms {C)D and 
{B)C are associated with the deterministic part. 

We then take the average over 6 and £ . Because of their Gaussian character, these 
averages are easy to carry out. Note that these averages give 



/exp J dt£ir,(t).(-0 Oi (t)\ = expt / dtdi! £ £ ^(t)-(6/o i (t)6/ 0j (t , ))-^(0 
and 

^exp J dtJ2 **♦(*)• [-90 I J\ ^e-^ fe2 ( t - t ')p T (k).^ (k, t')e^} ^ 

= expjEE J dtdt'ir t {ty ^ J dtdi' J J 6(t-t)6(t' _ t') e -^ 2 (*-f) e -^A 



(3.11) 



* 3 



X 



e ik.r lWe ik^ (i ')^( k; t)£J-(k', t')) 



(3-12) 



in (exp Lo)6>o£o an d contribute additional terms to the random part of the effective equation 
of motion, which is the term (B)A in Eq. ( |3.9| ). Bearing this in mind, we can easily write 
down the effective equation of motion. 
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C. Effective equation of motion and Markov approximation 



1. Effective equation of motion 

It is a simple matter to read off the effective equation of motion after the average over 
Langevin noises is taken. We can explicitly write down the expression for (B)A in Eq. ( p*.9| ) 
and attribute terms involving tensor to renormalized noise £ and the rest terms to renor- 
malized noise 0. (See Eqs. ( [3.15|) and Q3.16|) below.). Because the algebra is straightforward 
and slightly long, we will not write down the details here. The result is 



— 

dt ri[t) - Co^(t) 



1 dH 2B 2 2B 

5Wo —<Ji + Vi(t), 



Co 



(3.13) 



where 



r rt r N ~ 1 

Jo= e ik " r / dt'e^ k2 ^P T (ky € (k,0 - - E 

Jk J — CO L J0f) ._n 



95. 



°_ e -*Ti(tf) 



(3.14) 



is a vector field similar to the original solvent velocity field Eq. (|2.5b|) . The renormalized 
Langevin noises are still Gaussian, with variances 



OtifiOtf)) = (0 0i (t)0 0j (t')) + ^J 3 , 



(3.15) 



and 



^(k,t)^(k',t" 



(3.16) 



In these equations, for convenience, we have defined vectorial quantities J\ and J 2 and 
tensorial quantities J% and Their expressions are 

Ji= X - J dt'dld2(D(t,t', l,2))V i (t)E7(r 1 (t) - l)U(r t (t') - 2) 

+\Y,J dt'didld2(D(t',t,l,2))6(t - Or y (t,0-V i (t')t/'(r i (0 - l)£/(r,(f) - 2) (3.17) 



^ 2 = - / ^2(5^,^,1,2))^^)^^) - 1) 



+ \Y,J dt'dld2(B{t',t',l,2))e(t-t')r ij {t,t r )-V j {t')U{r j (f) - 1) (3.18) 
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J 3 = 5ij J dld2(B(t,t' 1 l,2))v i (t)V i (t')U(r i (t)-l)U(r i (t')-2) 

+2 J dmd2^B(i,t\l,2))6(t-r)T ij (t,t}-V j (t)V j (t')U(r J (t} - l)U(rj(lf) - 2) 
= 5 t] J dld2(B(t,t',l J 2))v t (t)V l (t')U(r t (t)-l)U(r t (t')-2) 

+2 J dtdld2(B(t,i, l,2))d(t' - fyViWTjitftfyViffiUiTi® - l)U(ri(t) - 2) (3.19) 

J" 4 =E / dU2(B(t,t',l,2))e- ik - T ^e- ik '- r ^ 

i 

xP r (k).V i (t)P T (k')-V i (t')^(r,(t) - 1)^(0 - 2) (3.20) 

(We have suppressed all the possible dependence of these J7"'s on time, monomer index, and 
so on.) 

Note that, although Eq. ( |3.14|) is essentially the same as Eq. fl2.5h| ), it would be incorrect 



to write a separate Navier-Stokes equation, similar to Eq. (|2.4cj) , for j7o, in which the full 
Hamiltonian is replaced by Ho, and claim that when the average over the quenched ran- 
domness is performed the effective equation of motion for solvent velocity is not changed. 
Recall that we have eliminated the solvent velocity field before the average is performed. In 
priciple, one can explicitly introduce an auxiliary field for the solvent velocity and perform 
the average over the quenched randomness to find effective equations of motion for both 
polymer conformation and solvent velocity. The effective equation of motion for solvent 
velocity field is presumably different from the original Navier-Stokes equation, though we 
expect that we would get the same equation as Eq. (|3.13|) once the velocity field is projected 



out. As we have mentioned in Sec. [0], our concern here is the equation of motion of the chain. 
Therefore, we will make no attempt on deriving an effective Oono-Freed kinetic equation for 
solvent velocity field. 

Following the conventional interpretation of MSR formalism |7|,|33|, time-correlation func- 



tions of dynamical variables {r^} and their auxiliary fields {f j} are response functions of the 
system to the perturbation associated with Causality requires that these response 

functions should be vanishing if any of the time variables of the auxiliary fields {? {\ is larger 
than all the time variables of the dynamical variables {r^}. Imposing this causality require- 
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ment upon Eqs. ( |3.10|) we can show that causality is also respected in the effective equation 
of motion, Eq. ( |3.13| ). 



2. Markov approximation 

The essential ingredients of the proof given in App. [A| are (1) elimination of solvent 
velocity field, and (2) Markov approximation. Physically, (1) means average over the solvent 
velocity and (2) means that the characteristic time scale for solvent dynamics is much shorter 
than the characteristic time scale for the motion of the chain. 

We adopt the same spirit of App. |A]. As the velocity field has been eliminated, we now 
apply the Markov approximation to further simplify the effective equation of motion. The 
Markov approximation is 

e -u k-(t-t') ^ 2 _ 
u k 2 

Within this approximation, the Oseen tensor 



T(r,r') = / J-P r (k)e ik -( r - r ') (3.22) 

appears, for example, 

IVM') = 2g 2 ( 8(t-t')T(r i (t),r j (t)) (3.23) 

Using Markov approximation we can largely simplify the effective equation of motion. 

Since we seek the result only to the lowest nontrivial order, that is, to the lowest order 
in B, go, and Bq, we can simplfy the expressions of J7"'s by dropping all the higher order 
terms. For convenience, we define 



p(R,t) = (X;*(r < (t)-R)) (3.24) 



G(t, t', 1, 2) = (£ - 1)^(0^(0 - 2)) (3.25) 
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C(t, t', 1, 2) = S(ri(t) - - 2)) (3.26) 

i 

where the bracket is the mean-field average defined by Eq. ( |B13|) in App. |B[ The mean-field 
kinetic equation, Eq. (|3.13|) , then reduces to 



i"W = -sS- s5 ? T(r " rj) -S + ' , ' (t) - (3 - 27) 

where the new Langevin noise n is Gaussian with variance 

^it)^{t')) = h(t - 0^1 + 2g%5(t - t)T(r t (t), R,-(t)) 



5 2 

-f 



2tfW dld2C(t ' t ' 1 ' 2) ^^) ar«(f) (3 ' 28) 

The effective Hamiltonian i? e // becomes 

# e// = H + B J2j dlp(r, t)U(Ti(t) - 1) 

"^E df / dlcBGfof , l,2)- ^ i( 2r 2) ^(t) ~ 1) (3.29) 

where t~ = t — e, e being an infinitesimal positive number. 



These equations have the form of the nonlinear generalized Langevin equation |44],[45 



However, the reason for the emergence of this nonlinear generalized Langevin equation is 



different from that in Refs. |4^j4q| . It is from the elimination of heterogeneity of the original 
system, rather than from the elimination of the degrees of freedom of surrounding "bath" 
molecules. Here the elimination of the surrounding "bath" molecules does not cause any 
memory effect because Markov approximation is used. 

In generalized hydrodynamics one can define frequency and wavelength dependent trans- 
port coefficients ||46|| . Similarly, in Eq. ( |3.29| ), the term containing heterogeneity B can be 
seen as giving a time-dependent excluded volume interaction. Alternatively, since this term 
is non-local in time, one can say that the heterogeneity introduces an excluded volume in- 
teraction along the world line, in addition to the usual excluded volume interaction B along 
the chain contour. 
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Inspecting these equations, we find that they are similar, albeit more complicated, to 
the effective kinetic equation obtained in spin glass dynamics |8]U. In these works a mem- 
ory function also appears in mean- field limit. In principle, one can expand the monomer- 
monomer interaction U in terms of dynamical variables r-j(t). This expansion will lead to 
a generalization of the p-spin-interaction spin-glass dynamics discussed in Ref. 0. How- 
ever, since one will get infinite p-spin-interaction terms, this approach does not seem to be 
practically feasible. 

We can now take the known result for homopolymer as a reference. When we set the 
chain heterogeneity B zero, Eq. (|3.27|) is expected to reduce to the same equation obtained 
from eliminating the solvent velocity field in the Oono- Freed kinetic equation (See App. 

| r ' W = "S-««STW*),r J W).g + rfM (3.30) 
fi(t)ff (<')) = l"5(« - OM + 2f$(t - OTMt), rj (t)) (3.31) 



where He is the Edwards Hamiltonian fl2.3|) . This equation reveals that, the lowest order 
approximation of p(R, t) must be the usual monomer density. Any higher order correction, 
after multiplied by B , must be first or higher order in B , gl/r/ , or B 2 /( . To the order of 
our calculation, these corrections can be dropped. Therefore, we can replace p(R, t) by its 
corresponding un-bracketed quantity: 

p(R,f) = X)(5fo(*)-R) (3-32) 

i 

Taking into account the symmetry of monomer indices, which should be restored when 
the self-consistent bracket is removed, the effective Hamiltonian ( |3.29| ) becomes 



H eff = H + ^J2U(r l (t)-r J (t)) 

z a 

"^E £" dt' J dld2G(t, t', 1, l)- 9U{ fc2') ^ U ( Ti ® ~ 1} (3 - 33) 
Since there is no corresponding known result for Bo = and B ^ (a heteropolymer 
chain in theta solution), we cannot find a lowest order approximation for G using the same 
method. 
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The quantity G(t, t', Ri, R2) in Eq. ( |3.33|) gives the response field of the chain at time 
t, when the chain conformation is S(ri(t) — Ri), to a perturbation, introduced by f ^ (t') at 
time t', when the chain conformation was 5(ri(t') — R2). (The perturbation can be explicitly 
introduced in the kinetic equation, in the same way as the perturbative magnetic field used 
in For simplicity it is not included here.) Causality again is respected and the integral 
of t' in the effective Hamiltonian is limited by current time t. 

The quantity C(t, t', Ri, R2) in the correlation function of the renormalized noise /Ltj, 
Eq. ( |3.28| ), gives the correlation function of chain conformations, between conformation 
S(ri(t) — Ri) at time t and 5(r,(t') — R2) at time t'. This is a natural dynamical generalization 
of the static order parameter proposed in Ref. J|]. 

Note that the effective Hamiltonian is a time-dependent Hamiltonian. The memory 
function G(t, t') makes the chain conformation at time t depend on all the past chain con- 
formations. This memory term makes the kinetic equation non-Markovian, and writing a 
Fokker-Planck equation for it impossible. It is well known that a non-Markovian process 
can be made Markovian by introducing new degrees of freedom (e.g. the exmaple in Ref. 
p7||). However, in our case, this would be equivalent to going back to the original Oono- 
Freed kinetic equations where the process was Markovian, because what has been done is 
exactly elimination of some degrees of freedom (heterogeneity and solvent) from the original 
Markovian process. Of course, a Fokker-Planck-like equation can be written for Eqs. (|3.27|) 



and (|3.33|) , but this would not give a closed equation. Instead, it gives an equation with 
infinite terms (e.g. a Kramers-Moyal expansion) ||31|| . 

On the other hand, a Fokker-Planck equation with memory effect has been developed by 



Zwanzig [48]. It is also shown in Ref. ||8] that, after certain approximations, the memory 
Fokker-Planck equation can be transformed into a generalized Langevin equation. This 
seems to allow us to transform the generalized Langevin equation obtained here back to a 
memory Fokker-Planck equation. However, it is not clear how to perform this backward 
transformation. 

These considerations make it clear that trying to write down a generalized (memory) 
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Kirkwood diffusion equation for the mean field equation of motion Eq. ( 3.27|) would be very 
difficult, if not impossible. Fortunately, although the relaxation spectrum is calculated in 
Ref. [|T7[] by perturbatively solving the Kirkwood diffusion equation, one can also derive the 
relaxation times directly from the kinetic equation. As it is clear from the above analysis, 
we should not try to find the corresponding generalized Kirkwood diffusion equation but 
rather derive the relaxation spectrum directly from the kinetic equation. This is done in the 
next section. 



IV. MODE RELAXATION SPECTRUM 



A. Transition to continuous model 



For convenience, we now rewrite the effective equation of motion in terms of continuous 



chain notation. We will use the standard notation |]T2|,|T9|: r being the contour length 
measured along the chain; c(r, t) the position of the point at r. A short distance cut-off 
a will be implicitly assumed in the following equations so that | r — r' \ > a for all double 
integrals over r and r' . The upper limit of contour length can be N — 1, according to Eq. 
( p.ip , or N, because N has been taken to be very large in mean field approximation. We 
will use N as the upper limit of contour length so that one can readily compare our result 
to the result of Ref. jl7| . 



The effective equation of motion in continuous chain notation is 



where 



M (r, *Mt', 0) = hs(t - t')5(r - t') + 2g 2 6(t - f )T(c(r, t) - c(r', If)) 

+ 2a 6{T - r) J dld2C{t ' t ' h2) ^^J) ^J) — (42) 

B 2 r , ft 



H e ff = H E - — 



[ dr /* dt' I dld2G(t, t', 1, 2y SU ( C } T / ,t \ s 2 V (c(t, t) - 1) (4.3) 

J J-oo J OC(tA') 
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HE = 2~J dT 1 ^ c(r ' t] |2 + i: I dTdr ' u « T i o - c ( r '' *)) ( 4 - 4 ) 

G(t, t', 1, 2) = /" dr(<y(c(r, *) - l)ic(r, 5 ( c (t'> *) ~ 2 )) ( 4 - 5 ) 

C(t, f, 1,2) = | dr(<y(c(r, *) - l)5(c(r', t) - 2)) (4.6) 

T(c(r,t) - c(r',t)) = [ -4p r (k)e !k,(c( ^ ) " c(T ' 1)) (4.7) 



From now on the potential U(c(r,t) — c(r', £)) will be taken as a hard-core repulsive 
interaction S(c(r, t) — c(r',t)) so that we can compare our result with Ref. ||17|| . For a 
general potential the corresponding equations can be easily written down. 

Eq. Q4.ip essentially describes a non-equilibrium process. In spin glass dynamics, as men- 
tioned above, a similar equation can be obtained through the same mean field approximation 
IHI. In Refs. [§,0, it is assumed that the process is stationary. While the newly developed 
out-of-equilibrium dynamics |05] is potentially able to address non-stationary dynamics, we 
will assume that in our case the process is stationary. This assumption gives 

G(t,t', 1,2) = G(t-t', 1,2) (4.8) 
C(t,t',l,2) = C(t-t',l,2) (4.9) 



B. Perturbation calculation 

The conventional method used to determine the spectrum of mode relaxation, especially 
the longest relaxation time, is to find the eigenvalues and eigenfunctions for the Fokker- 
Planck operator ||17|| . This approach is useful as long as a Fokker-Planck equation exists. 



However, as we have analyzed in the previous section, this is not the case here. 

Nevertheless, since what we want to find is the longest relaxation time, we can simply use 
perturbation method to find out how the kinetic coefficient is renormalized, with the same 
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spirit of Refs. [ 50|, 51]] . In principle, it is possible to develop a diagrammatic method, similar 
to the diagrammatics developed in Ref. ||52j| , to systematically calculate desired quantities. 
However, this would not make much sense because Eqs. ( |4.1| ) is valid only to the lowest 
nontrivial order. A straightforward perturbation calculation will be sufficient. 
Introducing Rouse coordinates $,(p,t), 

p-KT 



€(p,t) = J drc(r,t) cos 

and its Fourier transform 

€(p,u) = J dt£(p,t)e 
the equation of motion becomes 
G^Mtfau) = -^yf / drdr'dt jf e^ikcos 



N 



(4.10) 



iujt 



(4.11) 



~N~ 



x exp 



-k--£ I e-"«((p',J)Q r ,(T,T' 



piTT p'-KT 1 
COS — — COS ■ 



r] k 2 N 



N 



x exp 



B 2 



2C 2 V N 



k-E / e-^(/,^)g p »(r,r')lP T (k)^(p',^)e-^' 



pnr 
~N~ 



f dtdrdld2 C dill f e iujt G{t,t' ,1,2) -ik'ik cos 



xe- lk '- R2 exp 



xe ikRl exp 



N 



p' 



-iw't V 7rr 
cos 



N 



(4.12) 



and 



»(p,cu)»(p',cu')) = ~1 1 dtdre^+^cos^cos 



P TIT 



1^1 



N 


N 


piTT 


p'ttt' 






COS N 


C ° S N 


piTT 


p'ttt 






COS N 


C ° S N 



dld2C(t,t',l,2) 
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J6(c(r,t)-l)66(c(r,t')-2) 

6c(T,t) 5c(r,t') { ' ' 



where we have defined 



duj 
2^' 



(4.14) 



G-\ P ,u) = -iu + K p , K p = U^f, (4.15) 

Q p (r,r ) = cos— - cos^p. (4.16) 

Since we need to carry out calculations only to the lowest nontrivial order, we can iterate 
this equation once and drop all the higher order terms. The detailed calculation can be found 
in App. 0. It is shown in App. that the relaxation spectrum A p is renormalized by three 
terms at the lowest nontrivial order, 



SAP 



W^i^^^'^ (4 - l7) 



Ci 



rjo N N ' 
2QN 

Compared to the result obtained in Ref. ||17|| , Eqs. ( |4.17| ) and ( |4.18| ), renormalizations 



SA? = / dr | k Cos^cos^[G(r,0,k) - G(r, -u, k)] -k^, (4.19) 



due to excluded volume interaction and hydrodynamic interaction respectively, coincide with 
the results reported there. The additional renormalization due to the chain heterogeneity 
B is given by 5A p 3 \ It is found in Ref. |T7[ that the renormalizations from excluded volume 
and hydrodynamic effects are independent of each other, at the lowest nontrivial order. Here 
we see that the addition of heterogeneity preserves this independence: All three effects are 
not coupled. 

It has been shown that in spin glasses there exists a critical temperature below which 
the replica symmetric solution is not stable || . Below this transition temperature, there are 
infinite ground states so that the system displays non-ergodic behavior. This non-ergodic 
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behavior is responsible for many characteristics of the spin glasses below the transition 
temperature ||. It is proposed in Ref. || that, in an infinite system, the fluctuation- 
dissipation theorem (FDT) should be violated below the spin-glass transition temperature T c . 
Mathematically, the violation comes from the extra time-persistent term in the correlation 
and response functions. Physically, the breaking of FDT is caused by the non-ergodic 
behavior of the system. The response function contains not only the usual response obeying 
the FDT but also a time-persistent part corresponding to the crossing of the energy barriers 
between the ground states. Ref. || also shows that the breaking of FDT is needed in order 
to make the static mean-field solution of the Sherrington-Kirkpatrick (SK) model stable. 

It is clear that this dynamical consequence (i.e. breaking of FDT) is closely related to 
the static property (many ground states) of the SK model of spin glasses. Ref. [§J shows that 
a freezing transition can occur in heteropolymers at which only few conformations (states) 
dominate in the partition function. One might expect some dynamic consequences, and in 



fact it was suggested indeed ||T8|] that static freezing in heteropolymers is accompanied by 
that dynamic freezing. This assertion was made based on the phenomenological analysis 
which extends the REM to kinetics. The model used in (TTSJ employs an unphysical assump- 
tion that two states connected by one kinetic step have statistically independent energies; 
it is this assumption which leads to the conclusion about the glass transition in this system 
f53fl . Our microscopic analysis, Eq. fl4.19|) shows that, within the order of our calculations 
whether the FDT is broken or not, dynamically, there is no effect at all. This is because 
the G(t, 0, k) term exactly cancels the time-persistent part of the response field, if the tem- 
perature is below the transition temperature. Since the time-persistent part of response 
function is directly related to the degeneracy of the ground states, this result implies that 
this degeneracy is dynamically unobservable, at least to the lowest nontrivial order. This, 
in turn may mean that barriers between low-energy states are nonextensive (in N) in the 
studied model of heteroplymer. 

One may argue that the violation of FDT happens only when the system is infinitely 
large, N = oo, and, since only large N is imposed in the derivation of the effective equation 
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of motion, breaking of FDT does not happen and the above argument is not quite right. 
Nevertheless, what really matters is that the term G(r, 0, k) — G(t, —uj, k) makes all static 
effects at zero frequencies dynamically null, and corresponding static freezing phenomena 
may not have its have its dynamic counterpart in this system. 

Carefully checking the calculation that leads to the term G(t, 0, k) — G(r, —uj, k) in Eq. 
( |4.19| ), we see that it is rooted in the two exponential functions in the 5-dependent term in 
Eq. (|4.12|) . This makes it clear that the chain flexibility is responsible for the dynamically 
null result. Physically, this is a plausible result. The quenched randomness in Eq. Q2.1|) is 
quenched along the chain, but not in the space in which the heteropolymer is embedded. As 
we do not expect a dilute liquid mixture of magnetic and non-magnetic particles displaying 
any "frozen in" behavior observed in its solid mixture counterpart, e.g. spin glasses, we do 
not expect a full spin-glass-like behavior in heteropolymers because the randomness is not 
fully quenched in space. We probably can say that the correspondence between spin glasses 
and heteropolymers is statically more complete than dynamically. The analogy between 
spin glasses and biopolymers therefore has to be used with some care, especially in building 
ad hoc phenomenological models not supported by microscopic analysis. 



V. ON THE ANALOGY BETWEEN SPIN GLASS AND HETEROPOLYMER 

It is well known that there exists an elegant correspondence between magnetic system 
and self-avoiding random walk, e.g. a polymer chain with excluded volume effect. This 
correspondence can be shown perturbatively and analytically |]20| , |5^ , |55| . From the point 



of view of this complete correspondence, it seems that the analogy between spin glass and 
heteropolymer should be also complete. However, according the result obtain in the prvious 
section, one has to be careful when using this analogy. Therefore, it is necessary to think 
more about this analogy in the context of the correspondence. 

There are different ways to establish the correspondence between magnetic system and 
self-avoiding random walk [EDI. Here we will follow the analytical proof given in Ref. |5S . 
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Using the result obtained in Ref. |55[ we know that the model considered here corresponds to 
an Ising model with non-random nearest neighbor interaction and random four-spin coupling. 
More precisely, Eqs. ( |2.1| ) can be obtained from the n = limit of the following soft-spin 
Hamiltonian, 

n N N n R n 

# = EE J*™* + E + ^ (E O 2 ] (5.i) 

where r serves as a variable for Laplace transform (in getting the partition function for 
the walk) and its value needs not to be specified in our discussion. Further detail of this 
correspondence can be found in E) . 



Comparing this Hamiltonian with the one considered in spin glass dynamics |§ , in terms 
of soft-spin, 

H =\ E( r o^i _ - j - , ij) n ' n j + u E °f _ E h i a i ( 5 - 2 ) 

(ij) i i 

we see that the main difference is where the randomness enters. Eq. (|5.2|) has the quenched 
randomness being a two-spin coupling, while the random coupling Bij in Eq. (|5.1|) is in 
4-spin term. Note that the Hamiltonian (|5.1| ) is also different from the p-spin-interaction 
spin glass studied in Ref. , 

N 

H = - Jii-i p Vii---<7ip-J2 h i (J i ( 5 - 3 ) 

ii<i2'"<ip i=l 

(or, in terms of soft-spin, 

h = Eti7 a i + ua t) - P E Jii-ip^h ■ ■ ■ ^ P - P E h ^ ( 5 - 4 ) 

i 1 8=1 

). Although Eq. ( p.4| ) appears to be very similar to Eq. (|5.1|), the additional two-spin coupling 
Jij and the 0(n)-symmetry preserving term in Eq. ( p.l[ ) make the story entirely different. 
It is unclear, and seems unlikely, that the (qualitative) dynamic features obtained in Refs. 
|8||5[] have their counterparts in the dynamics of Eq. ( |5.1| ). In fact, we are not aware of any 
study of Eq. (|5.1| ) itself as a possible spin glass model. Perhaps Eq. ( |5.1|) being a spin glass 



model has not been studied, either statically or dynamically, before. 
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(Perhaps the 3 He- 4 He-aerogel system mentioned in Introduction is more relevant to the 
heteropolymer solution than spin glass system. We argued in the previous section that the 
randomness in heteropolymer solution is not completely quenched in space. In 3 He- 4 He- 
aerogel system, according to the picture proposed by Chan and coworkers 0, at low 4 He 
concentration, 4 He atoms can be seen as "partially quenched": Some are bound onto the 
aerogel, being completely quenched, while some are in freely moving superfluid phase. This 
interplay of quenched randomness and annealed randomness has been theoretically studied 
by Maritan et. al. ||56|| , using the so-called Blume- Emery- Griffiths (BEG) Hamiltonian [57 



Interestingly, the similarity between the BEG Hamiltonian [^,^1 an d Eq. ( |5.1|) seems higher 
than Eqs. Q and Q.) 

The above observation is only at the level of the Hamiltonian used in each model. In fact, 
more significant difference at the level of kinetic equation exists. In the Oono-Freed kinetic 
equation ( |2.4|) the dynamical variables are monomer positions {r^} and solvent velocity field 
u. On the other hand, the dynamical variables used in Refs. |S||| are spin components 
{<Tj}. When n — > 0, spins {<j,j} do not become monomer positions {r^} and solvent velocity 
field u. Although the effective equation of motion does share some common features, e.g., 
memory effect, renormalized noise, etc., with the ones found in spin glass study, these 
are the features common to the machinery, functional integral formalism plus mean field 
approximation. Having these common features still allows entirely different physics. As we 



argued at the end of Sec. |IV A| , our system is somewhat intermediate between the solid and 
molten phases of dilute magnets, and is perhaps more similar to 3 He- 4 He mixture in porous 
media than spin glasses. 

VI. CONCLUSION 



The main result of this paper is the effective equation of motion Eqs. ( |3.13| ) and ( |3.27| ) 
and the renormalized relaxation spectrum. 

We first ask how adaptive our model is. In polymer dynamics much work has been done 
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on homopolymers. There also have been some work on the dynamics of copolymers. As 
a first approximation, if the heterogeneity of the chain is not important, homopolymer is 
a reasonable model. It is seen in the above calculation that the lowest order effect of the 
heterogeneity parameter B is proportional to B 2 . This also indicates that this approximation 
makes sense. It is also known that the natural-occurred heteropolymers such as proteins are 
composed of at most 20 types of monomers (amino acids). In this sense the approximation 
of independent interaction energies B^ employed in this paper may be reasonable: It was 
shown 0| that the greater the number of monomer types is, the better the approximation 
of independent Bij is. 



Furthermore, in a recent simulation [p8[ , it is shown that in the early stage of the collapse 
of homopolymer, from coil state to globule state, some clusters form along the chain. One 
way to see this process is to view the chain, with clusters of various sizes, as an effective 
heteropolymer chain formed by monomers and clusters. This is plausible if one uses a coarse- 
grained viewpoint. One may argue that these clusters are not stable; they constantly form 
and annihilate. However, on average and in a coarse-grained sense, this constant birth- 
and-death of clusters probably can be ignored because the clusters seldom grow larger than 
certain scale, as the simulation shows. We may take the upper length scale of these clusters 
as a basic length scale in the model, and treat the clusters as heterogeneous monomers. 
These clusters may have different sizes and masses, which are also allowed in our model 
because sizes and masses of monomers do not appear in the Oono-Freed kinetic equation. 
In this sense, our model can also serve as a model for a homopolymer chain in poor solvent 
during its early stage of collapse. 



Ideally, similar to Ref. [fL7| , the next calculation would be applying renormalization 



group technique to find the scaling form of the relaxation spectrum. This does not seem 
promising, however. To perform this calculation, we need an ansatz for the w-dependence 
of G(lu, k), analogous to the one used in spin glass dynamics ||, as well as an ansatz for its 
k-dependence. Nevertheless, we do not expect that the renormalization group calculation, 
even if it is possible, will change the conclusion, regarding the dynamically null effect of the 
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time-persistent part of response function, reached in Sec. fVB , 

The most important development of the present model would be to incorporate chain 
compactisation, i.e. consider the dynamics of the heteropolymer globule. Another relevant 
and important aspect of polymer dynamics which may be studied along these lines is kinetics 
of coil-globule transitions pSJ which may be relevant also for protein folding. 
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APPENDIX A: DERIVATION OF KIRKWOOD DIFFUSION EQUATION 



We show here that the Kirkwood diffusion equation can be derived directly without the 



projection technique used in Refs. [12,16 



If the Hamiltonian is the Edwards Hamiltonian, the equation of motion derived from 
Eqs. fl2J) will be (cf. Eq. ( ^Wj ) 



—Ti(t) = --— -fr + g u${r i (t),t) + Qi (t) + SoulM*),*), 
at c,o OTi{t) 



<M)= / e 



ik-r 



y^ 1 dH E ik . r . (t /) 



Jk J — oo 



'k J — oo 

When Markov approximation is applicable 

e -u k 2 (t-t') _^ 

the equation of motion becomes 
dt ' 



Mt) = -c-odW)- 9 ~° 



At - 1') 



(Al) 



(A2) 



(A3) 



(A4) 



^ T(r,w ' r,(t)) 4 +0oiW+9o { e,k,r 'V pT(t) ^ (M) ' 

(A5) 
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Eq. flAjp is a Langevin equation with noise 8 (rj(t), t) = 0Qi(£) +fl , o u I( r i(^)^) satisfying 



e E {T i {t),t)6 E {v j {t'), 0) = 2Co5(t - + 2(7 2 T(r,(t), r,(t'))^(t - 



(A6) 



The Fokker-Planck equation corresponding to this Langevin equation is 

r d dH E 



+ 



n{r,},t) 



(A7) 



where P({rj},t) is the probability distribution function of chain conformation Eq. 
]) is exactly the Kirkwood diffusion equation. 



APPENDIX B: DERIVATION OF THE EFFECTIVE LAGRANGIAN 



Here we detail the calculation. 



Define 



Note that 



1. Ei#0?. term 



A(t, if, R, R') = E a i (t)a i (t')«y(r i (t) - R)5(r i (t') - R') 

i 

B(t, t', R, R) = E S(Ti(t) - R)5(r i (t') - R) 

i 

(7(t, R, R) = E ai(*)*(ri(t) - R)<Hr^') - R') 

i 

t', R, R) = E aitf)S(n{t) - R)S(Ti(lf) - R') 



A(t,t',R,R') = A(t',t,R',R) 
5(M',R,R) = £(t',i,R',R) 
C(t,t',R,R) = D(t',t,R',R) 



and they are all of order 0(N 2 ). Then, 



(Bl) 



(B2) 
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J2° 



- J dtdt'dld2d3dA x [A(t, t', 1, 3)B(t, t', 2, 4) + A(t, t', 2, 3)B(t, t', 1, 4) 
+A(t, t', 1, 4)77(t, t', 2, 3) + t', 2, 4)77(t, t', 1, 3) 
+C(t, t', 1, 3).D(t, 2, 4) + (7(t, 2, 3)L>(£, t', 1, 4) 
+C(t, t', 1, 4)£(t, 2, 3) + (7(t, 2, 4)L>(£, t', 1, 3)] 
xU(l-2)U(3-A) + 0(N) (B3) 

Let column matrix ip(t, t', 1, 2, 3, 4) and 16 x 16 matrix Si be 

fA(t,t', 1,3)\ ( B{t,t',2,A)\ 



^ A {t,t', 1,2,3,4) = 



^ c (t,f, 1,2, 3, 4) = 



A(M',2,3) 
A(*,t',l,4) 
U(M',2,4)/ 

/C'(t,t',l,3)\ 
C(M',2,3) 
C(t,t',l,4) 



Mt, i 7 , 1,2,3, 4) = 



^(t,f, 1,2, 3, 4) = 



5(t,t',l,4) 
B(t,t', 2,3) 
V B(t,t', 1,3) J 

( D(t,t',2, 4) \ 

£(M',M) 
D(t,t', 2,3) 

\D(t,t', 1,3) J 



ip(t,t', 1,2,3, 4) 



Si 



(B6) 



/^ A (M',1,2,3,4)\ 
^(t,t',l,2,3,4) 
^(M',1,2,3,4) 
VVd(M',1,2,3,4)7 
(where 1 4 is a 4 x 4 unit matrix), then 

]T C?2 = y dtdt'dld2d3dAi) T (t, t', 1, 2, 3, 4)SW(*, t', 1, 2, 3, 4)17(1 - 2)17(3 - 4) + O(N) 

(B7) 



( 1 4 \ 

1 4 

o u 

V o o i 4 o ) 



(B4) 



(B5) 



2 - Ej^j Cij term 

The calculation for Ci? can be performed in a similar way. However, because of the 
symmetry of C and D, there is an ambiguity in expressing ®%j as a sum of products of 
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C and B, or of C and B, or a mixed type. This should not make the final result different, 
as we expect the same symmetry can be used to convert C and D. 
If we express @ij in terms of products of C and B, then 



°ij = 7 / dtdt'dld2d3d4il) T {t, t', 1, 2, 3, 4)S 2 i/)(t, t', 1, 2, 3, 4)Z7(1 - 2)<f(* - t') + O(iV) 



(B8) 



5 2 = 



(B9) 



/0 0\ 

o u o 

1 4 
\0 0/ 

where the overbarred delta function 5 is used to keep track the fact that the integral over 
time variables, hence, the 5-function, has to be done after the integral over space variables 
has been performed. 



3. Gaussian transform and mean-field approximation 



Combining the results from the previous two sections, (e La ) b in Eq. ( |3.8|) can be written 



as 



exp{ / dtdt'dld2d3d4[— ^ T S^U{1 - 2)U{3 - 4)-A/> T S 2 ipU{l - 2)6{t - £')] + 0{N)\ 
J &Co 4£ J 

(BIO) 

Introducing variables Qi(l, 2, 3, 4) • • • Qi6(l, 2, 3, 4) and performing a Gaussian transform, 
the generating functional Z 9o ^ becomes 

ZBoto = J {VTiHViiHVQi} exp[~^Q T AQ + Q T ip + L + O(N)] (Bll) 



where 



16 



(B12) 



4=1 
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Q is the column matrix formed by Q\ ■ ■ ■ Qie and A is the coefficient of the quadratic term 
in Eq. flETop . 

Using steepest decent method the mean-field approximation gives 



Q°(M',1,2,3,4 
Q°(t,t', 1,2, 3, 4 
Q°(f,*',l,2,3,4 
Ql(t,t', 1,2,3,4 
Q° 5 (t,t', 1,2,3,4 
Q° 6 (t,t', 1,2,3,4 
Q° 7 (t,t', 1,2,3,4 
Q° 8 (t,t', 1,2,3,4 
Q° 9 (t,t', 1,2,3,4 
Q° l0 (t,t', 1,2,3,4 
Q° n (t,t', 1,2,3,4 
Q° 12 (t,t', 1,2, 3, 4 
Q° l3 (t,t', 1,2,3,4 
Q° u (t,t', 1,2,3,4 
Q? 5 (M',1,2,3,4 
Q%(t,t', 1,2,3,4 



B 2 
B 2 



—^U(l-2)U(3-4)(B(t,t',2,4) 
4 so 



4Co 2 
5 2 

4Co 2 
B 2 



U{l-2)U{3-4)(B{t,t',l,4) 



772^(1 - 2)17(3 - 4)(S(t, t', 2, 3) 



U(l - 2)U(3 - 4){B(t,t',l,3) 



4Co 2 

-^[/(l-2)f/(3-4)(A(t,t',l,3) 

-^[/(l-2)[/(3-4)(A(t,t',2,3) 
4 s,o 

4Co 2 

4Co 2 

^U(l-2)U(3-4)(D(t,t>,2,4) 

= -^f/(l-2)f/(3-4)( J D(t,t',l,4 
£? 2 



U(l -2)U(3-4)(A(t,t',l,4) 



-j^Uil - 2)U{3 - 4)(A{t, t', 2, 4) 



4Co 2 
B 2 



U{l-2)U{3-4)(D{t,t',2,3 
U(l-2)U(3-4)(D(t,t',l,3 



4Co 2 

^U(l-2)U(3-4){C(t,t',l,3 

^U(l-2)U(3-4)(C(t,t',2,3 
4 so 

f- I U(l-2)U(3-4)(C(t,t>,l,4 
B 2 



U{l-2)U{3-4)(C{t,U,2,4 





2)5(t - 


t')(C(t,t', 1,3)) 




2)5(t - 


t'){C{t,t', 2,3)) 




2)(!(t - 


t')(C(t,t', 1,4)) 




2)5(t - 


t')(C(t,t',2,4)) 




2)5 (t - 


t')(B(t,t',2,4)) 




- 2)5{t - 


-t')(B(t,t', 1,4)) 




- 2)5{t - 


-t')(B(t,t', 2,3)) 




- 2)5{t - 


-t'){B(t,t', 1,3)) 



where the angular bracket means 
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/{Pr. t }{Pr t }(---)e^e^ 

The mean-field solution for Q® are determined self-consistently from these equations. 
Therefore, the effective Lagrangian L e = Lq + (Q°) t i/j is 

B 2 r 

L e = L +-j / dtdt'dld2d3dAU{l - 2)U(3 - 4)[(A(t, t', 1, 3))£(t, t 1 , 2, 4) 
Co ^ 

+(B(t, t', 2, 4)>A(t, f , 1, 3) + 2(C(t, t\ 1, 3))D(t, t', 2, 4)] 

2fin 



+— - / tfttft / dld2d3d4<y(t-t / )?7(l-2)[((7(t,t / ,l,3))S(t,*',2,4) 
Co J 

+{B(t,t',2,A))C(t,t',l,3)}, (B14) 



APPENDIX C: FIRST ORDER CORRECTION OF RELAXATION SPECTRUM 



The equation of motion in Sec. [IV B| has the following form 

G \p, w)£(p, u) = flip, co) + F[& p, uj] (CI) 
Substitute the zeroth order solution 

C% 1 (p,u)tf a >(p,u>) = i*(p,u) (C2) 
into it, the first order solution is 



^\p,co) = G (p,u)»(p,u) + G (p,u)F[& \p,u] (C3) 

This solution contains terms to the lowest nontrivial order. Since £ enters in .F through 
exponential function, further iteration does not change the result to the lowest nontrivial 
order. Therefore, to the order we seek for, the solution is 

£(p,w) = Go(p,u>)n(p,u>) + Go(p, uj)F[^°\p, u\ + higher orders in (C4) 

To find the relaxation spectrum, we find the average (over /x of the following tensor 
product 
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£(p, u)£(p', u/)= G (p, uj)G (p', u/)h(j>, uj)fi{p', u') 

+Go{ P ,co)G ip',u;')[^i P ,u;)F^ip',u;') + ^p',u')F^(p,u;)} (C5) 

where F^{p,u) = .F[£^,p, u;], and the term containing two F^(p,u) has been dropped 
because F^(p,u) is at least the lowest nontrivial order. 
Let 



B n (t, t') 



(C6) 



We find 



£(°>(p, ^)exp[z J|k- £ | # e--^(^, W ')^(r, 
v p' w 

G (p^)Go(p,-^)|^ke^g p (r,r')e- fc2B -^') 

C (0) (p, «0exp[<^k. £ e--"^(°)(p", o;")Q P «(r, /)] 

2 

G (p,uj)G (p, -u) — [2n5(uj + u')6 PtP >l - 
Co 



(C7) 



(C8) 



Co iV 



(p, c)exp[zy^k. £ £ £(p', cos ^(e-' 



ia/i g— iu> t 



Gob, w)G (p, -o;)^y pk{e™ - e**) cos 



As noted in Sec. [IV A| we confine ourselves to the stationary case. Therefore, 



(CIO) 



G(t,t',k) = G(t-t',k) 



(Gil) 



C(t,t',k) =C(t-t',k) 



(C12) 
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where, assuming G(t, t', 1,2) = G(t, f, 1 - 2) and C(t, t', 1, 2) = C(t, t', 1-2), i.e. transla- 
tional invariance, 



G(t,t',k) = yrfRG(t,t',R)e^ k - R (C13) 

C(M',k) = / rfRC(t,t , ,R)e- lk - R (C14) 

Defining 

G(r,w,k) = j dte luJt G(t,k)e- k2B ^ (C15) 

(B*(T,t,t') = B*(T,t-t')) and 

C(r,o;,k) = J dte^C^e'^ 3 ^ (C16) 

then (<i is the dimensionality) 

(£(p,o;)€(pV)> 



27T(5(W + uj')G (p, uj)G (p, -lu) 

-Co P ' p 9o NJ dTdT L r^ 2 C ° S iV C ° S N 6 1 

5 2 2 f f p-KT p'W ~, , . 

+ 2^ivi ^kkcos — cos— C(r,u;,k) 

+~G (p,u;)ll(p,p', u) + ~G (p', -u)K(p',p, -u) (C17) 

iV (, iV (,o 



^iJdr jf t[G(r,0,k) -G(r,o;,k)].kkkcos^cos^ (C18) 

Note that this equation is obtained after higher order terms are dropped. Therefore, it does 
not mean that c<j) is a Gaussian process. 
Define 



5A« = "f | / drdr' / k ^ cos ^Q p (r, /)e- fe2 ^(-') (C19) 
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5I& = -^l(^) 2 fdrdr' [ i^cos^cos^e-^^') (C20) 
p rjo N N J h k 2 N N K ' 

8Af\u) = J dr / k Cos^cos^z[G(r,0,k) - G(r, c, k)]-k^ (C21) 

The scalar contraction of the tensor product gives (We take p = p' mode only because of 
orthogonality of modes.), 

(£(p s w)-£(p,u/)) = G (p,u)G Q (p, -u)(fi(p,uj)-fi(p,uj')) 

x[l - G (^)(5A« + 5Af + Mf(u;)) 

-G (p, -^)(5A« + 5Af + 5Af(-^))] (C22) 

where (/Lt(p, u)-fji(p, ct/))o is the value obtained from Eq. Q4.13| ) by using zeroth order solution 
Eq. 

Therefore, 

[G^{p,u) + 5A p (cu)][G 1 (p,-co) +5A p (-u)](£{p,un(p,u')) = (,x(^MpV)>o (C23) 

5A p (uj) = 5A« + 6A® + 6A® {u) (C24) 

The relaxation spectrum is therefore renormalized 

A p -> A p + 5Ap(cj) (C25) 

We note that the first two corrections, SA p ^ and 5A^' due to excluded volume effect and 
hydrodynamic effect, respectively, are the same as those found in Ref. [fTTjl . 
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